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Motivated by its role as a central pillar of current theories of dynamics of spin ice in and out 
of equilibrium, we study the single-ion dynamics of the magnetic rare earth ions in their local 
environments, subject to the effective fields set up by the magnetic moments they interact with. 
This effective field has a transverse component with respect to the local easy-axis of the crystal 
electric held, which can induce quantum tunnelling. We go beyond the projective spin-1/2 picture 
and use instead the full crystal-held Hamiltonian. We hnd that the Kramers vs non-Kramers nature, 
as well as the symmetries of the crystal-held Hamiltonian, result in different perturbative behaviour 
at small helds (< 1 T), with transverse held effects being more pronounced in Ho 2 Ti 207 than in 
Dy 2 Ti 207 . Remarkably, the energy splitting range we hnd is consistent with time scales extracted 
from experiments. We also present a study of the static magnetic response which highlights the 
anisotropy of the system in the form of an off-diagonal g tensor and we investigate the effects of 
thermal huctuatlons in the temperature regime of relevance to experiments. We show that there is 
a narrow yet accessible window of experimental parameters where the anisotropic response can be 
observed. 


I. INTRODUCTION 

The properties and behaviour of spin ice materials, 
Ho 2 Ti 207 (HTO) and Dy 2 Ti 207 (DTO) among oth¬ 
ers, are deeply rooted in the characteristic single-ion 
anisotropy of rare earth (RE) magnetism. The two low¬ 
est energy states (degenerate at single-ion level) are sep¬ 
arated by a large energy gap (> 200 K) from the other 
excited states, thns projecting the system onto an effec¬ 
tive spin-1/2 space at low temperatures. The lowest en¬ 
ergy donblet has moreover a strong easy-axis anisotropy, 
which is responsible for its classical Ising-like behavionr^ 
(for a recent detailed discnssion, see also Ref. 2). These 
properties jnstify modelling the magnetic moments as 
classical Ising spins with a local easy axis. The rich ther¬ 
modynamic behaviour of spin ice systems can be largely 
accounted for by the physics of the ground state dou¬ 
blet combined with the pyrochlore lattice structnre and 
exchange and dipolar interactions: frnstration leads to 
an extensively degenerate gronnd state^, topological or¬ 
der, and an emergent gange symmetry hosting magnetic 
monopole excitations^’^. 

This thermodynamic model of spin ice was later pro¬ 
moted to a dynamical one by introducing an experi¬ 
mentally inspired®^® single spin-flip time scale® (see also 
Ref. 10). This choice was motivated by the experimen¬ 
tal observation of a well-defined microscopic time scale in 
the magnetic response of these materials, which appears 
to be largely temperature independent in the regime of 
interest. Such dynamical modelling of spin ice proved 
reasonably snccessfnl at capturing the experimental re¬ 
sponse and relaxation properties, and triggered a new 
research direction into the behaviour of these systems 
out of equilibrium - an interesting and highly tuneable 


setting that combines topological properties, kinematic 
constraints, emergent point-like quasiparticles and long- 
range Conlombic interactions^^’^®. 

A temperatnre-independent microscopic spin-flip time 
scale is typically associated with quantnm tnnnelling un¬ 
der an energy barrier that the ion has to traverse in order 
to reverse its magnetic polarisation. Understanding this 
behavionr clearly requires that we go beyond the single¬ 
ion ground-state doublet (spin-1/2) approximation, and 
we investigate the role of possible qnantum perturbations 
that may be responsible for the tnnnelling dynamics. To 
date, such understanding appears to be lacking in the 
literature. 

The work presented here is a step at gaining insight 
into the quantum single-ion dynamics in spin ice HTO 
and DTO. Specifically, we focus on spin-spin interactions 
as a sonrce of quantnm fluctuations. The exchange and 
dipolar fields acting on a given ion due to others in the 
system have both a longitudinal and a transverse com¬ 
ponent with respect to the local easy axis. The latter 
acts as a transverse field in the effective Ising model. 
We study in detail the effects of snch transverse field on 
the single-ion behaviour, obtaining the resnlting energy 
splitting (namely, inverse characteristic time scales) and 
anisotropic response, both at zero as well as finite tem- 
peratnre. 

There exists a concrete motivation for stndying the 
specific case of an exclnsively transverse magnetic field, 
a setting which at first sight seems to reqnire flne-tnning 
the longitudinal component to vanish. This may not ap¬ 
pear straightforward in spin ice, a dense assembly of large 
rare earth moments interacting via a long-range and ge¬ 
ometrically complex dipolar interactions, Eq. (1). 

However, spin ice is no stranger to snch fine-tuning. 
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It is now well understood how the geometry of the py- 
rochlore lattice conspires with that of the dipolar inter¬ 
action to ensure that the longitudinal total field on each 
spin is, to a good approximation, equal for all spins in all 
ground states^^, which are exponentially numerous^^ and 
in general not related by any symmetry transformation. 

Similarly, a pointlike defect in a spin ice ground state, 
known as magnetic monopole^ has an energy independent 
of its location, as long as it is spatially well separated 
from other monopoles. As the spatial displacement of 
a monopole proceeds via the flip of a single spin, this 
spin must be subject to a vanishing longitudinal field - 
otherwise the excitation energy in the system, encoded 
in that of the monopole, would change as the monopole 
moves. 

Therefore, our study can be thought of as providing a 
picture of the quantum mechanics underpinning the mo¬ 
tion of an isolated monopole defect in a ground state of 
spin ice. The properties of such mobile monopoles are 
currently subject to both experimental^®^^^ and theoret¬ 
ical work^®. 

For the purpose of the present paper, we approxi¬ 
mate the exchange interactions by their classical form. 
Namely, we consider the interaction Hamiltonian 


H = 


-jJ2s,-s, 

(ij) 


(.ij) L 


Si ■ Sj 3(5'i • rij){Sj ■ rij) 


...13 


( 1 ) 


where i,j label the sites of the pyrochlore lattice; Si = 
Ui'Li, with fji = ±1 and are the (4 inequivalent) unit 
vectors pointing from one tetrahedral sublattice to the 
other; J and D are the exchange and dipolar coupling 
constants, respectively; is the nearest neighbour dis¬ 
tance on the pyrochlore lattice; and is the distance 
between the two sites i and j. Within the approxima¬ 
tion of this Hamiltonian, the action of all other ions on 
a given one is an effective magnetic field, 

H = Y,Bes{^)-S, 

i 

Bes{i) = — J Sj 

3, (ij) 

+ Drl^ ^ 
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whose strength and direction were studied in Ref. 19. 
Here we thus limit ourselves to considering the action of 
an applied field on the full single-ion Hamiltonian, be¬ 
yond the customary projection to its lowest lying states. 

We And that the Kramers vs non-Kramers nature of 
HTO vs DTO results in a different perturbative be¬ 
haviour at small fields, whereby in HTO the ground 
state doublet splits at second order in the applied held 


S, 


3(5, 


strength, whereas DTO only splits at third order, as illus¬ 
trated in Fig. 5. One therefore expects transverse fields in 
HTO to be more effective at inducing quantum tunnelling 
dynamics than in DTO (at small fields, < 0.5 T). Using 
degenerate perturbation theory, we provide an analytical 
understanding of this difference in behaviour in terms of 
symmetries of the OF Hamiltonian. Remarkably, the en¬ 
ergy splitting range we And is consistent with quantum 
tunnelling time scales observed in experiments^. 

We also present a detailed study of the static mag¬ 
netic response to a transverse field, which highlights 
the anisotropy of the system. We And interesting reso¬ 
nances as a function of the in-plane direction of the field, 
where off-diagonal components of the g-tensor become 
non-vanishing (namely, a purely transverse held in the 
cry-plane induces a longitudinal response along the local 
z axis). 

We investigate the effects of thermal fluctuations in 
the temperature regime of relevance to experiments. We 
And that in thermal equilibrium much of the anisotropic 
response averages out up to rather large fields a few 
Tesla) for temperatures as low as a hundred milliKelvin. 
Nonetheless, signatures of the anisotropic response in 
spin ice could be experimentally observed at low tem¬ 
peratures in fields ~ 10 T. Our results further support 
the robustness of the classical easy-axis Ising approxima¬ 
tion for the single-ion behaviour in spin ice, while at the 
same time helping to quantify its limit of validity. 

We stress that all of the above features can only be 
grasped via a full description of the single-ion Hamilto¬ 
nian capturing the complexity of its interaction with the 
other spins / environment. They cannot be understood 
(but at most added at an effective level) if we limit our 
modelling projectively to the lowest CF levels. 

The paper is organised as follows. Sec H introduces 
the full single-ion crystal-field (CF) Hamiltonian for rare 
earth ions Ho®+ and Dy®+ in spin ice. Sec. HI investi¬ 
gates the effects of a magnetic field at zero temperature, 
with specific focus on fields transverse to the local easy 
axis. We use exact diagonalisation (Sec. HI A) as well as 
degenerate perturbation theory in the limit of small fields 
(Sec. HIB). Thanks to the large CF energy scales typical 
of these systems, perturbation theory is indeed valid well 
into the range of held strengths of experimental interest. 
Finally, Sec. IV discusses thermal effects in the relevant 
temperature range and studies the behaviour of the re¬ 
sulting single-ion magnetic susceptibility, and in Sec. V 
we summarise and discuss our results. 


II. CRYSTAL-FIELD OF SPIN ICE RE®+ IONS 

The general formula for spin ice pyrochlore oxides 
is A 2 ^B 2 ~*’Oy“, where the A and B species are rare 
earth (RE) and transition metal (TM) cations, respec- 
tively^°^^^. The structure is given by the space group 
BdSm featuring two sub-lattices that interpenetrate each 
other and consist of networks of corner-sharing tetrahe- 
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dra. In HTO and DTO the A magnetic sites host, re¬ 
spectively, the Ho^+ and Dy^+ ions, while the B sites are 
occupied by non-magnetic Ti^"*' ions. 

The local point group symmetry for the RE^+ ions 
in magnetic pyrochlore oxides is a trigonal (see 

App. A). This is schematically shown in Fig. 1 and it 
accounts for the arrangement of the eight oxygen ions 
(yellow spheres) surrounding the rare earth ion (green 
sphere). The oxygen sites are distinguished in two main 
subclasses according to their position with respect to the 
central RE-site: the 01 sites and the 02 sites. The strong 
axial alignment of the 01 ions (above and below the cen¬ 
tral RE^+ ion) drives the classical Ising-like anisotropy 
typical of spin ice materials. The 02 ions, displaced in 
equilateral triangles lying in parallel planes transverse to 
the easy axis of the 01 ions, are responsible for the an- 
tiprismatic character of the symmetry. 



Figure 1. The crystal-field environment of a ion in 

a magnetic pyrochlore oxide. Ho®+ is used for concreteness. 
Panels (a), (b) and (c) show respectively a tilted, side and top 
view of the same structure. For clarity, panel (a) addition¬ 
ally displays the six surrounfing Ho^+ ions and their second 
axial oxygens. The Ti'*"'" ions are arranged in an hexagon 
coplanar with the ion in the centre. The edges of the 

triangles connect the coplanar 02 oxygens: three above and 
three below the plane of the RE®"*" and Ti'^^ ions. The two 
green planes shown are parallel to each other and contain the 
respective triangles of 02 ions. The antiprismatic arrange¬ 
ment of the six 02 gives the Dsd point-group symmetry. The 
two remaining Of ions, aligned along the (111) axis with the 
central RE®"*", drive the local Ising anisotropy. 


The crystal-field Hamiltonian of a rare earth ion in a 


D^d symmetry can be conveniently expressed as^^’^^ 


HcF ^B^oOl + B^d^ + Bldl 
+Bld^o + Bldl + Bldl 


( 2 ) 


where the Stevens operators together with the re¬ 
spective parameters B^ determine the CF spectrum and 
eigenfunctions of each compound. Following the general 
convention^^, the are such that the k = 0 opera¬ 
tors are g-polynomials of only diagonal operators 3^, Jz, 
while those with k > 0 include also fc-powers of the lad¬ 
der operators J+, . A list of the matrix elements of the 

Stevens operators in the |J, Mj) basis, where J,Mj are 
the quantum numbers for the total angular momentum 
and its projection along the local (111) axis respectively, 
is given in Ref. 25, and can be straightforwardly obtained 
from their operator expressions in App. B. The crystal- 
field parameters for HTO and DTO are listed in Table I. 



HTO (meV) 

DTO (meV) 

Bo 

-7.6 X 10"^ 

-1.6 X 10“" 

Bo" 

-1.1 X 10"® 

-2.3 X 10“® 

B| 

8.2 X 10“® 

1.6 X 10"® 

Bo® 

-7.0 X 10"® 

6.5 X 10"® 

Bg® 

-1.0 X 10“" 

9.9 X 10"® 

Be® 

-1.3 X 10"" 

1.0 X 10"" 


Table 1. The crystal-field parameters (in meV) for the Hamil¬ 
tonian in Eq. (2) obtained from Refs. 26 and 27 (see also 
App. C). 


The CF Hamiltonian can be diagonalised to obtain the 
CF states. The spectrum is, in general, made of multi- 
plets and singlets since the Stark splitting, induced by 
the crystalline electric fields, removes only partially the 
2J + 1 degeneracy of the ground state multiplet. The 
spectrum of HTO (Fig. 2a) features five singlets and six 
doublets, while the spectrum of DTO (Fig. 2b) is only 
made of doublets. This discrepancy is due to Kramers 
theorem forbidding singlets in spectra of atoms with an 
odd number of electrons (Ho^+ has n = 10 electrons in 
the 4-/ shell, while Dy^+ has n = 9). The order of mag¬ 
nitude for the energies, however, is roughly the same, and 
the ground state is a doublet in both. The energy gap 
between the ground state doublet energy and the first 
excited level is in excess of 200 K. 

Two possible basis eigenfunctions for the ground state 
doublets are displayed in Fig. 3, showing that they can 
be well approximated by the fully polarised states IV'o) ~ 
\Mj = J), I'f/ji) « \Mj = —J). This illustrates the strong 
anisotropy along the local quantisation axis in both sys¬ 
tems. 
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Figure 2. crystal-field spectra for HTO (a) and DTO (b), 
respectively. The spectrum of HTO features both doublets 
(solid lines) and singlets (dashed dotted lines). In meV, bot¬ 
tom to top, the series of doublets is: 0, 21.96, 25.99, 59.59, 
71.51, 76.80 while the one of singlets is: 20.42, 27.71, 69.36, 
69.94, 80.52. In contrast, since Dy'^'*’ is a Kramers ion, DTO 
features only doublets. These are eight in total: 0, 25.23, 38.0, 
38.21, 51.75, 77.49, 87.65, 89.16. Note the thicker line just be¬ 
low 40 meV is not a quadruplet, but rather it corresponds to 
the two doublets 38.0, 38.21. 


III. EFFECT OF A MAGNETIC FIELD 





(b) DTO 


The degeneracy of the crystal-field spectra is removed 
in the presence of a magnetic field B: 

= ?icF — j ■ B. (3) 

In this equation, /re = eli/2me is the Bohr magneton 
(e and mg are, respectively, the charge and the mass of 
the electron) and gj is the Lande factor for the RE^+ ion 
with total angular momentum J {gj = 5/4 and gj = 4/3, 
respectively, for Ho^+ and Dy^+). 

In the following, we use the local coordinate system 

xo = ;^(l,l,-2), yo = ^(-1,1,0), zo = ^(1,1,1), 

( 4 ) 

with respect to the global axes X, Y, Z of the cubic py- 
rochlore unit cell (see Fig. 4), with the Zq axis conve¬ 
niently pointing along the high-symmetry direction of the 
crystal-field Hamiltonian^®’^^. 


A. Exact diagonalisation 


Figure 3. A possible choice of basis wave functions, |^/)o) 
(solid lines) and (dotted grey lines), for the ground state 
doublet of Ho®'*' (a) and Dy^+ (b). The wave functions have 
been obtained by diagonalising the crystal-field Hamiltonian 
in Eq. (2) with the crystal-field parameters given in Table I. 
Insets: the same shown over a narrower vertical range. 



A longitudinal field along the local easy axis leads to 
conventional Zeeman splitting linear in field strength and 
selects of one of the two polarised states in Fig. 3. At sim¬ 
ilar field strengths, this is the field direction that results 
in the largest energy splitting due to the anisotropy. 


Figure 4. The local coordinate frame xo,yo,zo (red arrows) 
used to describe the transverse magnetic field B (green arrow) 
and its direction angle (f) in Eq. (5). The atoms are shown in 
the same top-view as in Fig. Ic. 
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If the longitudinal component vanishes and a purely 
transverse field component is present, then the chosen 
polarised basis states split into symmetric “bonding” and 
antisymmetric “anti-bonding” combinations. A polariza¬ 
tion in the plane perpendicular to the easy axis is however 
opposed by the anisotropy and this competition results 
in unusual effects that will be discussed in the following. 

The coupling of the total angular momentum to the 
transverse magnetic field can be written in terms of lad¬ 
der operators as 

J • Bx = i |B^| (e-*^J+ + e*^J_) . (5) 

In the local coordinate system, </> is the angle of the field 
with respect to xq in the plane transverse to the easy axis 
zo (see Fig. 4). 



10^ T). This illustrates the strength of the energy scales 
set by the crystal-field and provides a reference for mag¬ 
netic field values that can be considered a small pertur¬ 
bation. 

At lower fields, when the two competing terms in 
Eq. (3) have comparable energies, the response of the 
system becomes anisotropic. This anisotropy is much 
stronger for Ho^+ than for Dy^+ and, for (p = 30° -I- n 60° 
with n integer, it leads to resonances (due to level cross¬ 
ing between Eq and Ei) shown in Fig. 5 (red dotted- 
dashed line) and in Fig. 6a. 
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Figure 6. Ground state splitting for both HTO (a) and DTO 
(b), as a function of the angle p at four particular values 
of the transverse field (from bottom to top in each figure): 
B = 0.03 T, B = 0.55 T, B = 11.31 T, and B = 228.33 T. 
For HTO, the two strongest field values resonate and make 
the splitting close aX (p = 30° -I- n60° (n integer). 


Figure 5. Splitting of the ground state doublet under the in¬ 
fluence of a purely transverse magnetic field. The red curves 
correspond to the non-Kramers behaviour of HTO, while the 
blue curves correspond to the Kramers behaviour of DTO. 
Note the y axis is dimensionless to allow a consistent com¬ 
parison of the two systems. The different curves correspond 
to fixed angles of the transverse field: 0 = 0° (solid curve), 
0 = 10°, (short-dashed), 0 = 20° (dotted) and 0 = 30° 
(dotted-dashed); angles of 0-|-n 120°, with n integer, give ex¬ 
actly the same curves because of the OF trigonal symmetry 
of the 02 ions. The two long-dashed straight lines show the 
limiting behaviours at very high fields - Larmor precession, 
Eq. (6) - and at low fields for HTO - degenerate perturbation 
theory, Eqs. (7a) and (8). 

The dependence of the splitting of the ground state 
doublet AEqi vs the magnitude of the transverse field 
is shown in Fig. 5 for both HTO and DTO. For very 
large fields the anisotropic effect of the OF environment 
becomes negligible and the magnetic moments undergo 
simple Larmor precession with frequency wp given by 

AAoi = feuL = gjHB |B| . (6) 

Due to the strong crystal-fields in HTO and DTO, 
such regime is clearly experimentally unattainable {B > 


Finally, at fields of the order of 1 T or less, the ion 
enters a perturbative regime where AEqi is given by the 


following power-laws: 



AAoi = |B|' 

for HTO, 

(7a) 

AAoi= 44(0) |B|° 

for DTO. 

(7b) 


Ho3+ is not a Kramers ion and features some singlets 
in its unperturbed energy spectrum. These are responsi¬ 
ble for the quadratic behaviour (low field asymptotics in 
Fig. 5), whose coefficient 

4^)0 = 2.68x10-°^ ( 8 ) 

can be obtained analytically from perturbation theory 
(see Sec. HIB). Dy°+ instead is a Kramers ion and all 
unperturbed energy levels are doublets. As we explain in 
the next section, this causes the quadratic correction to 
vanish identically, leading to a cubic dependence on the 
applied field. Fitting the corresponding asymptotic low- 
field behaviour in Fig. 5, we obtain the angle-dependent 
coefficient 

«DTO (<^) = 6.8 X 10"'^ (1 -f A cos (60)) (9) 
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with A = 0.114. 


B. Perturbation theory 

We can gain insight into the low field behaviour by 
using (degenerate) perturbation theory on 

n=Ho-XV, ( 10 ) 

where Ho = Hcf is the CF Hamiltonian in Eq. (2) and 
the perturbation V = £cf J • B/|B| corresponds to the 
Zeeman energy in Eq. (5), tuned by the dimensionless 
parameter A = (/j/tb |B| /£cfi where fcF is an arbitrary 
reference energy scale, e.g., related to the CF bandwidth. 
It is useful to introduce as the (unperturbed) CF 

eigenstates with energy En ^ (n = 0, ...,2J). 

The splitting of the RE^+ ground state doublet is given 

by 


AF;oi = AV(Eo.o-V^i.i)"+4|Eo.i| 


+ A' 


\ 




\fe>i 


AK 


(0) 

O/c 


E 


Vo.fcI4,i 

(oT 


k>l 


( 11 ) 


up to second order in A. In this expression Vn,m = 

{'ip^n^\V and a4°) = 4°) - 

Firstly, we notice that both HTO and DTO have Vo,o = 
Vi,i = Vb,i = 0 (see App. E), and therefore the first order 
contribution vanishes identically. 

In order to evaluate the second order contribution, we 
need to consider matrix elements of the transverse field 
perturbation V between the ground state doublet and the 
excited states. The symmetries of these matrix elements 
reflect the symmetries of the crystal-field environment. 

In App. E we discuss the different contributions in de¬ 
tail. We find two different behaviours for the excited 
states that form doublets. Some of them (type A) have 
identically vanishing matrix elements with the ground 
state doublet, Vb,n = W.n = 0 for n belonging to type 
A, and they trivially do not contribute to the splitting 
at second order. The other doublets (type B) have non¬ 
vanishing matrix elements which satisfy the following re¬ 
lations: 

|R0,m|" = iRl.m+ll" (12) 

= \Vo.m+l\^ (13) 

Vo,mVm,l + 4o,m-|-lKn+l,l = 0, (14) 

where ) and |^/>Ei) are the two eigenstates belong¬ 
ing to an excited doublet of type B (see Eq. (E6) and 
Eq. (ElO) in App. E). These relations imply that also 
doublets of type B do not contribute to the splitting 
at second order. (Details of the A- and B-type doublet 
wavefunctions and their matrix elements with the trans¬ 
verse field operator are given in Tables IV-V in App. E.) 


These results hold for both DTO and HTO. The former 
only has doublets (of either type A or B) in the spectrum 
due to Kramers degeneracy and no splitting occurs at 
second order. Indeed, the second order term in Eq. (11) 
generally reduces to the sum of the contributions from 
the singlets alone (see Eq. (E7) in App. E). Fig. 5 clearly 
shows that a non-vanishing third order contribution does 
exist, which we extract by fitting, Eq. (9). 

HTO, on the contrary, has some singlets amongst its 
excited states, which give a non-vanishing second-order 
contribution to the splitting. This can be readily com¬ 
puted, Eq. (8), and is in excellent agreement with the 
slope found from the numerical simulations (see the cor¬ 
responding long-dashed straight line in Fig. 5). We notice 
that the third order contribution has an angular depen¬ 
dence on (j) which is absent at second order. 

It is interesting to notice in Fig. 5 that the cubic power- 
law found for DTO persists up to ^ 100 T. In contrast, 
the quadratic power law in HTO begins to break at fields 
of the order of 0.1 T, depending on the in-plane angle </>, 
holding up to almost 10 T for (p = 30°; for all other 
angles, the cubic term becomes clearly dominant in the 
range from 1 T to 10 T, with an angular dependence 
similar to the one for DTO. 


C. Doublet splitting and time scales 

Let us compare the observed ground state doublet 
splitting in Fig. 5 with experimental magnetic relaxation 
time scales in spin ice^. The latter are typically of the 
order of 1 ms (at least in DTO), which corresponds to an 
approximate energy splitting of 10“^ — 10“® K. 

In order to estimate the former, one needs typical 
values for the exchange and dipolar transverse field 
strength. Ref. 19 suggests the range 0.1 — IT. Using 
Eqs. (7a), (7b), (8) and (9), we find that this corresponds 
to splittings in the range of 10“® K to 10“® K. 

This rough theoretical estimate is consistent with 
the experimental value. Whilst further investigation is 
clearly needed, the result is nonetheless suggestive that 
internal fields generated by exchange and dipolar inter¬ 
actions can in principle be responsible for (single-ion) 
quantum spin-flip dynamics in spin ice. 


D. Anisotropic response to a transverse field 

The strong single-ion anisotropy plays a crucial role 
also in the magnetostatic behaviour of the RE^+ ions at 
zero temperature. This is illustrated in Fig. 7, panels a-f, 
as a function of the angle (p and of the field strength | B |, 
where (Jq.) = {iplJalip), and a = x,y,z label the three 
components for the local coordinate system in Fig. 4. 

Both HTO and DTO acquire negligibly small values of 
(Ja;) and {Jy) for fields up to 10 T (see also Sec. IVA),. 
Moreover, we observe a sizeable (zero temperature) re¬ 
sponse in (Jz) to a purely transverse field, signalling non 
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Figure 7. Panels (a-f): Expectation values for the three components of the total angular momentum J in the ground state 
of the Hamiltonian in Eq. (3) with purely transverse magnetic field: {Jo) = {^| JcIV'), ol = x,y,z, as function of the angle (j) 
and the strength |B| of the field in logarithmic scale. For both HTO (left) and DTO (right), the x,y components are negligible 
for fields up to 10 T. In contrast, the rriz components feature a sizeable periodic dependence on the angle (p below 10 T. This 
is a manifestation of the strong axial anisotropy characterising the ground state of the spin ice RE^+ ions. Note the different 
response in the two systems: DTO features a smooth angular dependence which becomes asymptotically constant in the low 
field limit (from 10 T down to the lowest fields), whereas the oscillatory behaviour in HTO is more abrupt and its amplitude 
descreases from the saturated value reached at approximately 10 T, down to zero at low fields. Panels (g-h): Finite temperature 
behaviour of the expectation value (J^) = Tr(Jap)/Tr(p) at T = 0.5 K for HTO (left) and DTO (right). The Boltzmann weights 
from the density operator average the two (lowest-energy) states with opposite polarisation along zq. 

























vanishing off-diagonal components of the g-tensor. Of 
course, for high enough fields (|B| ^ 10 T), all expecta¬ 
tion values tend to the angular dependence of the Larmor 
regime, as expected when the Zeeman energy dominates 
over the CF Hamiltonian. 

The main difference between HTO and DTO is in 
the behaviour of (J^) below 10 T. In HTO, for fields 
IT < |B| < lOT, {Jz) oscillates rather abruptly with 
respect to the angle <j) between the saturated values —8 
and 8; for fields below 1 T the amplitude of oscillation de¬ 
creases and it becomes vanishingly small at low fields. On 
the contrary, in DTO the angular dependence is smoother 
and approaches constant (maximum) amplitude for fields 
below 10 Tesla; the amplitude however never reaches sat¬ 
uration. The period of oscillations is the same in both 
HTO and DTO (120 °), but we observe a phase difference 
of 60°. 

We notice that the behaviour of (J^) in Fig. 7 is con¬ 
sistent with the arrangement of the oxygens surrounding 
the rare earth ions. Depending on the (j) angle in the 
xo,yo plane (Fig. 4), the magnetic field can take three 
inequivalent high symmetry directions: either towards an 
oxygen that lies above the plane (0 °-|-n 120 °), or towards 
an oxygen that lies below the plane (60 ° +n 120 °), or else 
precisely in between two oxygens (30 ° -I- n 60 °), where in 
all cases n is an integer. The first two directions corre¬ 
spond to the maxima and minima of = gjHB {Jz)-, 
respectively. The latter direction corresponds to nodes 
where ruz vanishes. Curiously, as we noted before, the 
sign of ruz in the first two cases switches between DTO 
and HTO. We notice that this switching is highly depen¬ 
dent on the precise values of the CF parameters used. 


IV. FINITE TEMPERATURES 


panels a-d), we focus our discussion on {Jz)- Its be¬ 
haviour as a function of (j) and |B| is shown in Fig. 7, 
panels g-h, for T = 0.5 K. 

We find that the anisotropic response survives at inter¬ 
mediate fields, in between a high and a low field thresh¬ 
old. The high field threshold is the (temperature in¬ 
dependent) onset of Larmor precession. The low field 
threshold instead is set by the ground state doublet split¬ 
ting (Fig. 5). The low field threshold is temperature 
dependent, namely ~ for HTO and ~ (i.e., 

more easily observed) for DTO, according to the results 
in Sec. HI. 


B. Magnetic susceptibility 

The susceptibility of the a-component of the magnetic 
moment with respect to the /3-component of the applied 
field B is given by 


d{L) 


(15) 


At high temperatures, we expect the system to behave 
as an ordinary paramagnet, whose zero-field magnetic 
susceptibility (per spin) is given by the Curie law 

X% = , ( 16 ) 

where g? = J( J -I- 1). Therefore, it is convenient to 
define the dimensionless quantity 


Xa/9 _ SfeeT d{Ja) 
X^ 9j1^bJ{J+1) OB 13 


The splitting between ground and first excited state 
can be very small at low fields (see Fig. 5), far smaller 
than any temperature of experimental interest. Since 
the two states originate (adiabatically) from the split¬ 
ting of the ground state CF doublet, they will, in general, 
preserve the (symmetric) property of being polarised in 
opposite directions (that is, the equivalent behaviour of 
Fig. 7, panels e-f, for the first excited state - not shown - 
is simply opposite in sign with respect to that for the 
ground state). In the absence of a longitudinal field, 
thermal averages between the two will therefore cancel 
out the single-ion moment. 

A. Magnetic moment 

At finite temperature T, (Jq) = Tr( jQ/5)/Tr(/5), where 
p = jg density operator in the microcanoni- 

cal ensemble, "H is the Hamiltonian in Eq. (3), and fee is 
the Boltzmann constant. 

Since (J^,) and {Jy) take on negligible values at applied 
fields below the (trivial) Larmor threshold (see Fig. 7, 


whose behaviour is shown for /3 = a in Fig. 8. All curves 
exhibit Curie behaviour at (unphysically) high temper¬ 
atures. (Only Xxx and Xzz are shown, as Xyy behaves 
analogously to Xxx-) 

The behaviour of Xzz at low temperatures is per¬ 
haps most remarkable. It exhibits a Curie-like interme¬ 
diate temperature regime where « 2.5 x^i for 

both HTO and DTO. At high temperatures, this regime 
crosses over to the expected Curie law at an (approx¬ 
imately) field-independent threshold T ^ 10^ K, set by 
the CF energy gap between the ground state doublet and 
higher excited states. Below this threshold, the system is 
effectively projected onto its ground state doublet. The 
magnetic response is thus enhanced since these two states 
carry the largest magnetic moments of all CF levels. 

In presence of a finite applied magnetic field, as is the 
case in Fig. 8, one trivially expects a lower threshold to 
the effective spin-1/2 Curie behaviour when the temper¬ 
ature becomes smaller than the (linear) Zeeman split¬ 
ting between the two levels. The system then crosses 
over to a regime where the susceptibility is temperature- 
independent, but finite, corresponding to a small residual 
polarizability in the ground state. 



9 



(a) HTO (b) DTO 

Figure 8. Logarithmic plots of Xaajx^, with a = x, z, as a. function of temperature T and in presence of a static applied 
field. When the temperature is lowered the xa;-component deviates from the Curie law by decreasing approximately linearly 
while, the z 2 -component exhibits an intermediate (higher) plateau ~ 2.5). (Note that Xaa/x^ oc T corresponds, 

following Eq. (17), to Xaa being T-independent. The dashed line in each panel illustrates a linear behaviour as a guide to the 
eye.) Each component is shown for three values of applied fields: 0.001 T, 0.01 T, and 0.1 T. For x^^^ the three curves overlap 
almost perfectly, signalling that the susceptibility is field-independent below 0.1 T. For Xzz, the three curves overlap only for 
sufficiently large (field-dependent) threshold temperatures; in the main text we discuss how this behaviour is directly related to 
the different ground state splittings opened in the crystal-field spectrum by the applied fields. The insets show (Jx) vs field Bx 
at different temperatures, demonstrating a linear regime up to at least 1 T. The temperature-independent susceptibility below 
10 K is reflected in the perfect overlap of the magnetisation curves (Jx) vs Bx at these temperatures {Xxx ~ C /Sks, with 
(7 = 0.015for HTO and C = 0.02DTO). 


Interestingly, Xxx displays a similar behaviour, in spite 
of the fact that the splitting between the two lowest- 
lying states in a transverse field is now much smaller than 
temperature. In this case the temperature-independent 
regime extends all the way to T ^ 10 K. 

V. SUMMARY AND DISCUSSION 

We have presented a detailed study of the single-ion 
behaviour in spin ice HTO and DTO in presence of an 
applied magnetic field, based on the full description of 
the single-ion crystal-field Hamiltonian. We have con¬ 
sidered both zero and finite temperature, and focused in 
particular on the case of a field transverse to the local 
easy axis. We find that the Kramers vs non-Kramers na¬ 
ture of HTO vs DTO results in a different perturbative 
behaviour at small fields. 

We also present a detailed study of the static mag¬ 
netic response to a transverse field, which highlights the 
anisotropy of the system. We find that as a function of 
the in-plane direction of the field, off-diagonal compo¬ 
nents of the g-tensor become non-vanishing (namely, a 
purely transverse field in the xy-plane induces a longitu¬ 
dinal response along the local z axis). 

Within the classical exchange and dipolar Hamiltonian 
approximation, the action of all other ions on a given one 
is an effective magnetic field, whose strength and direc¬ 


tion were studied in Ref. 19. The transverse component 
of this effective field can be thought of as a potential 
source of quantum dynamics in spin systems, and the cor¬ 
responding ground state splitting studied in the present 
paper corresponds in this view to an inverse character¬ 
istic time scale. It is then remarkable to notice that - 
in spite of these simplifying assumptions - the resulting 
time scales for HTO and DTO are consistent with the 
ones observed in experiments '. 

There are a number of natural directions for future 
work. One is towards a yet more microscopic picture, go¬ 
ing beyond the effective field approximation we have em¬ 
ployed here, by determining the actual exchange Hamil¬ 
tonian for example via a superexchange calculation. An¬ 
other lies in considering the interplay of the spin degrees 
of freedom in the same spirit as we have considered the 
coupling to an external field. For the case of magnetoe¬ 
lastic couplings, a simple calculation in this spirit was 
reported in Ref. 28. 

Indeed, the issue of coupling to non-magnetic degrees 
of freedom is of relevance given the remarkably long 
millisecond-timescale of the spin flip and the small split¬ 
ting of the ground doublet. These are below 10~^K in 
temperature units, well below the scale at which experi¬ 
ments are conducted. Understanding the spin tunnelling 
process in the presence of a coupling to the ‘hot’ environ¬ 
ment is therefore an interesting exercise, close in spirit 
to the study of molecular magnets, which may be of rel- 















10 


evance to some of the unexplained features of the slow 
low-temperature dynamics of spin ice. 

Many of our results can be tested experimentally, per¬ 
haps best in heavily Y-diluted systems in an externally 
applied magnetic field, to reduce the added complexity of 
spin-spin interactions. Fig. 7 shows, that the anisotropic 
response of a RE ion in spin ice could be observed at 
temperatures ~ 100 mK under externally-applied fields 
~ 10 T. At lower fields, the induced magnetic moment is 
much lower, but, as the insets of Fig. 8 show, it could still 
be detectable, for example by muon spin rotation. Such 
experiments could provide a quantitative validation of 
the present description, which is a crucial step towards 
gaining further insight into the quantum dynamics of spin 
ice materials. 
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Appendix A: Oxygen environment 

The crystal-field interactions, i.e., the Stark effect due 
to the negative charges of the oxygen ions, deeply affect 
the single-ion quantum states. The physics dictated by 
the crystalline fields of the oxygens is so fundamental 
that, in the context of magnetic pyrochlore oxides, often 
it is preferable to use the expression A 2 B 2 O 6 O', instead 
of A 2 B 2 O 7 , simply to emphasise the role played by the 
oxygens according to their crystallographic and ligand 
character. Referring to a given RE^+ (A) site, e.g., a 
ion, the oxygens are arranged around it in an anti- 
prismatic fashion which is often referred to as a distorted 
cube (see Fig. 1 in the main text). 

The level of distortion is, however, huge compared to 
an ideal cube, as the two 01 oxygens form a linear O-A-0 
stick oriented normal to the average plane of the remain¬ 
ing six 02 oxygens arranged in triangles above and below 
the central A ion. The A-01 and A-02 bond distances 


are different: the former, ~ 2.2 A, is amongst the shortest 
bonds ever found in nature; the latter can vary depending 
on the compound, although in general it is between 2.4 
and 2.5 This implies that each RE^+ ion is char¬ 

acterised by a very pronounced axial symmetry along the 
local ( 111 ) axis which joins the two centres (01 sites) of 
the tetrahedra, through the magnetic ion sitting at the 
shared vertex (see Fig. la). The axial symmetry is af¬ 
fected by the anti-prismatic arrangement of the 02 ions 
with respect to the central RE^+ ion. These, as shown in 
Fig. 1 in the main text, are grouped in triangles lying on 
planes, above and below the RE^+ ion, which are parallel 
to each other. 


Appendix B: Stevens operators 


The Stevens operators in Eq. (2), expressed in 
terms of the angular momentum operators, can be writ¬ 
ten as: 


dl=ZJl-3^ 

64 = 35 J4 25 _ 30J2j2 - 6 J 2 -k 3J^ 

6g = 231J® -k (735 - 315J2) Jt 

/ . , N . (Bl) 

-k f 294 - 525J2 -k 105J^ j 

- 60J^ -k 40J^ - 5J® 

Ol=l{ {Jl + ii) , (11 J," + 59 J, - 33^1) } ; 

6l = l[jl + Jl) , 

where J± = Jx ^iJy and the anticommutator |A, = 

Ab + BA. 

Since each 6^ is a function of J = (^Jx,Jy,J^, the 
total angular momentum operator of the magnetic ion, 
the single-ion CF states can be conveniently expressed in 
terms of | J, Mj), where J, Mj are the quantum numbers 
for, respectively, the total angular momentum and its 
projection along the local (111) axis. A list of the matrix 
elements of the Stevens operators in the | J, Mj) basis is 
given in Ref. 25. 


Appendix C: Derivation of the Stevens crystal-field 
parameters 

The interaction between a magnetic RE^+ ion and its 
surrounding crystalline environment is usually described 
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starting from the simple Hamiltonian 


= - ^ |ei|VcF(f*), (Cl) 

i 

where Vcf represents the crystal-field potential from to 
the surrounding ions acting on the electrons in the un¬ 
filled shells of the central RE. Each i-th electron feels 
a potential VcF(ri) = VcF{ri,di,(j)i) at position ri. To 
study the crystal-field interaction it is convenient to make 
use of spherical coordinates centred on the RE site be¬ 
cause of the spherical symmetries of the electrons of an 
atomic system^®. It is customary to write the Hamilto¬ 
nian in terms of tensor operators. The tensor operator 
for the i-th. electron is 

= (C2) 

and obeys the same transformation rules as the spherical 
harmonics. In terms of these, the CF Hamiltonian for a 
magnetic ion in a crystalline D^d point-group symmetry 
reads^^ 

ncF = Blcl + + Bt{Ct - CU) 

+ Blct + Rf(Cf - (763) + BliCt + C \), 

(C3) 

Here the sum over the 4-/ electrons (X]r=i) omitted 
together with the index i for simplicity. 

The B^ parameters encapsulate the effect of the sur¬ 
rounding charges. Eq. (C3) is thus an alternative no¬ 
tation to the one based on Stevens operators, Eq. (2). 
The latter is convenient for RE^+ ions, where \J,Mj) is 
a good basis for the quantum states of the correlated 4-/ 
electrons, because they are explicit functions of the an¬ 
gular momentum operators J^, J^, J_|_, J_^^. The B^ are 
related to the B^ by means of the following expressions: 




for g > 0 . 


(C4) 


The are the factors outside the square brackets [... ] 
in the list of tesseral harmonics in Cartesian coordinates 
in Table IV of Ref. [25]. The 6 k (with k = 2,4,6; 62 = 
otj^d^ = /3j, 06 = 7 j) calculated by Stevens for different 
RE ions^® are given in Table VI of the same Ref. [25]. 
In Table H we reproduce the values for Oj, /3j, yj for the 
two magnetic ions Ho®+ and Dy®®" in spin ice materials. 

Experimental techniques based on inelastic neutron 
scattering are the most suitable to measure accurately 
the crystal-field energies in real compounds. From these 
measurements a reliable estimate of the CF parameters 
can be inferred beyond the level of accuracy allowed by 
the point-charge approximation®®’®^’®®. 

The crystal-field energies and parameters common in 
the literature of spin ice materials are based mainly on 



Ho®+ 

Dy®+ 

Oij 

-1 

-2 

450 

315 

Pj 

-1 

-8 

30030 

135135 

IJ 

-5 

4 

3864861 

3864861 


Table II. The 9k values (respectively a,j, fij, yj for A: = 2, 4,6) 
for Holmium and Dysprosium trivalent ions®®. 


the experiment presented by Rosenkranz et al. in Ref. 26. 
There, the neutron scattering measurement of all the CF 
energy levels allowed a complete parametrisation of the 
Hamiltonian in Eq. (C3). The full list of the B^ parame¬ 
ters for HTO found in that reference is reproduced in the 
first column of Table HI. Similarly for DTO, the second 
column of Table HI gives the suggested in Ref. 27 
as an interpolation of the values known for Ho 2 Ti 207 
and Tb 2 Ti 207 . However, to the best of our knowledge 
no neutron scattering experiments have been carried out 
successfully to determine the CF parameters of DTO. 
The corresponding R^ parameters, listed in Table I for 
the Hamiltonian Eq. (2) in the main text, follow from 
Eq. (C4) and Table HI. 



HTO (meV) 

DTO (meV) 

Ro 

68.2 

51.1 

Ro 

274.8 

306.2 

to 

83.7 

90.5 

Ro 

86.8 

100.4 

Rf 

-62.5 

-74.4 

R| 

101.6 

102.9 


Table III. Crystal-field parameters for the tensor opera¬ 
tors formalism. The parameters for HTO have been measured 
by means of inelastic neutron scattering in Ref. 26. The ones 
for DTO were derived as an interpolation of the parameters 
known for Ho 2 Ti 207 and Tb 2 Ti 207 in Ref. 27. 


Appendix D: Degenerate Perturbation Theory 

The results in Sec. HI A, namely Eqs. (7a-ll), follow 
from conventional degenerate perturbation theory (see 
e.g.. Ref. 31). In this Appendix we outline the main 
steps of the derivation for convenience. This will help 
understand the role that the symmetries of the unper¬ 
turbed CF Hamiltonian play in determining the pertur¬ 
bative behaviour, as discussed in the main text and in 
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detail in App. E. (The interested reader can find all the 
details of the calculations in Ref. 32.) 

We use the notation 

H = Ho + \V, (Dl) 

where A is a small (real) parameter tuning the strength 
of the perturbation V, and where the eigenstates and 
eigenvalues of Hq are 

Ho IV'f) . (D2) 


is not the case for singlets, which are present in the spec¬ 
trum of non-Kramers HTO, where splitting occurs at sec¬ 
ond order and Eq. (D4) is in good agreement with the nu¬ 
merical solution at sufficiently small values of the pertur¬ 
bation parameter A (see main text). On the other hand, 
Kramers theorem forbids the appearance of singlets in 
DTO resulting in a vanishing second order splitting. 

Appendix E: Matrix elements of the perturbation 


Here Hq = Hcf, in Eq.(2) of the main text, and V is 
the applied magnetic field, see Eq. (5) and Eq. (10). We 
focus on the case of interest where the first two energy 
levels are exactly degenerate Of course, in 

this case the choice of basis for the ground state doublet, 
|'0o°^) and is not unique. 

Expanding both eigenstates and eigenvalues of the per¬ 
turbed Hamiltonian in powers of the parameter A, one 
obtains the contributions to the GS doublet splitting or¬ 
der by order. The form of the contribution at a given 
order depends of course on whether the two levels did or 
did not split at lower order. 

For notational convenience, it is useful to define the 
matrix elements 


1. First order 

The first order contribution is bookwork, 

^(Eo^o-Vi4)2+4|KoaP. (D3) 

However, (main text and App. E) both DTO and HTO 
have Vb.i = Vo.o = Vi,i = 0 so that the degeneracy is not 
resolved at first order in A. 


2. Second order 

We focus on the case of interest where Vq.o = = 

= Vb,i = 0. After a few lines of algebra one obtains 
that the second order (GS doublet) contribution to the 
splitting takes the form 



|Eo,fcP- |Ei,fcP 


AE 


(0) 

k 


2 

+ 4 


E 

fe>i 


Vo,kVk,i 

AE^°'> 


2 


(D4) 


where AE^^'^ is the energy difference between the GS dou¬ 
blet and the fc*** excited state of the unperturbed Hamil¬ 
tonian Hq. The summation over A: > 1 spans all unper¬ 
turbed states other than the GS doublet. 

As discussed in the main text and in App. E, we find 
that all excited states that form doublets in the unper¬ 
turbed spectrum amount to a vanishing contribution to 
the second order GS doublet splitting in Eq. (D4). This 


In this Appendix we give details of the calculation of 

the matrix elements V where \ip^) 

are the eigenstates of the unperturbed GF Hamiltonian, 
and the dimensionless operator 

V = J+ -k 6+*“^ J_ (El) 

represents the applied magnetic field (perturbation) 
purely transverse to the local quantisation axis of the 
RE-ion (see Eq. (5) in the main text). Namely, the oper¬ 
ator V = £cf J ■ B/|B| in Sec. HI A relates to Eq. (El) 

via V = 2Vj£cF- The perturbative regime corresponds 
to field values within the initial power-law behaviour of 
the splitting observed in Fig. 5. 

The quantum states for HTO and DTO are represented 
in Table IV and Table V, respectively. Each state |^) is 
given by a superposition |'0) = |Afj). Reading 

from the left, the first two columns account for and 

li’A+i)- These are called A-doublets of the GF spectra 
to underline their different structure compared to \'4 ’b'^) 
and IV’s+i)) the H-doublets listed in the fifth and sixth 

columns. The ground states, |'0 q°^) , belong to the 

A-type doublets (explicit values of the coefficients are 
given in the captions of each table). 

In the third and fourth columns, the states 

and V \'ipA+i)j obtained by applying the perturbation V 
to the A-doublets, are given to facilitate the calculation 
of Vb,m and Vi,m, i-O., the coupling between the ground 
state doublet and the other GF states. The perturbed 
states in the third and fourth column are expressed in 
terms of coefficients defined as: 

V \Mj) = jtt, \Mj + 1) + j-^ \Mj - 1), 

(E2) 

The only depend on the angle (j) of the field in 
Eq. (El) and on the quantum numbers J,Mj. Further¬ 
more, from the general properties of the ladder operators, 
we have 


Jmj = j-iMj±i) ’ (E3) 

which leads to characteristic symmetries of the Vi,m el¬ 
ements that are key to determine the behaviour of the 
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Table IV. The coefficients for the decomposition of the single-ion states of HTO with respect to the angular momentum eigen¬ 
states \M,j). A generic state at the top of a column of coefficients Cmj is given by a superposition |i/)) = l^j)- 

It is useful to distinguish between A-type and B-type doublet states, as well as singlets s and s'. The ground state doublet 
IV’o°^) 1 is of the A-type; the coefficients are {ag = 0.982,05 = 0.156,02 = 0.065, o_i = 0.071, o _4 = 0.049, o _7 = 0.006} 

and account for the strong Ising anisotropy characteristic of the spin ice single-ion physics. The second and third column from 

the left are the first order terms obtained by applying V to an A-type doublet. 


leading orders of the perturbative splitting of the ground 
state doublet. Since HTO is a non-Kramers system, Ta¬ 
ble IV also shows two kinds of singlets in the last two 
columns on the right. 


1. HTO 

The crystal-field spectrum of HTO is made of 5 singlets 
and 6 doublets (see Fig. 2a). As summarised in Table IV, 
there are two types of doublets, A and B, and two types 
of singlets, s and s'. 


(V'i+ll V = 0-708 (^J-7 - Js ^ + 


a. A-doublets 

The ground state doublet in HTO, |V’o°^) ^^d is 

made up of A-type states, as defined in Table IV. 

It is then immediate to prove that Vo,o = 0 and 

Vi 1 = 0 since in general IK = 0 and 

(V'a+iI ^ IV'a+i) = 0- Namely, the first and the third 
column (and the second and fourth) of Table IV have 
trivially vanishing overlap. 

On the contrary the first and fourth (second and third) 
columns have a priori non-vanising overlap. In order to 
see that again Vi^o = hop = 0 one ought to consider the 
explicit form of the matrix elements 

-405 (^jZ4 -j5^+ 0-102 (^jZi - , (E4) 


which vanishes because all elements within round brack¬ 


ets cancel out, according to Eq. (E3). 
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DTO 

IV’™) 

IV’aIi) 

V|v>f) 


IV’™) 

IV’sil) 


15/2) 

C15/2 

0 - 15/2 

0 

0 

0 

0 


13/2) 

0 

0 

•/15/2 015/2 

/15/2 0 - 15/2 

& 13/2 

6 - 13/2 


11/2) 

0 

0 

j^/2 09/2 

~^g/2 0-9/2 

611/2 

—6-11/2 


9/2) 

C9/2 

— 0 - 9/2 

0 

0 

0 

0 


7/2) 

0 

0 

jg /2 O9/2 

~^g/2 0 - 9/2 

67/2 

— 6 _ 7/2 


5/2) 

0 

0 

js/g 03/2 

^3/2 0 - 3/2 

65/2 

6 - 5/2 


3/2) 

C3/2 

0 - 3/2 

0 

0 

0 

0 


1/2) 

0 

0 

h/2 “3/2 

h/2 “-3/2 

61/2 

6-1/2 


-1/2) 

0 

0 

J- 3/2 «-3/2 

7 - 3/2 “ “ 3/2 

6-1/2 

—61/2 


-3/2) 

c-3/2 

“03/2 

0 

0 

0 

0 


-5/2) 

0 

0 

•/-3/2 0 - 3/2 

-7I3/2O3/2 

6 - 5/2 

—65/2 


-7/2) 

0 

0 

j-g/2 0 - 9/2 

jtg /2 “9/2 

6 - 7/2 

67/2 


-9/2) 

a-9/2 

09/2 

0 

0 

0 

0 


-11/2) 

0 

0 

jIg/2 0-9/2 

jZg /2 “9/2 

6-11/2 

611/2 


-13/2) 

0 

0 

•/- 15/2 0 - 15/2 

“ 7 - 15/2 O15/2 

6 - 13/2 

—613/2 


-15/2) 

C- 15/2 

—015/2 

0 

0 

0 

0 


Table V. The coefficients for the decomposition of the single-ion states of DTO with respect to the angular momentum 
eigenstates \Mj). The CF states for DTO are all doublets. These are either of type A: IV’a+i) 

second column, respectively; or of type B: IV’s^) £^nd in the fifth and sixth column, respectively. The ground state 

doublet belongs to the A-type states; the coefficients are { 015/2 = 0 . 983 , 09/2 = — 0 . 171 , 03/2 = 0.044, a_ 3/2 = 

0.044, o_ 9/2 = 0.008, o_i 5/2 = 0}. The second and third column are the first order terms obtained by applying F to a pair of 
A-doublet states. 


This shows not only that Fo,o = ^ 1,1 = hi,o = ^ 0,1 = 
0 , accounting for the vanishing first order splitting in 
Eq. (11) in the main text, but also that all matrix el¬ 
ements coupling the ground state with any other A- 
doublet of the CF spectrum have to be null. Summaris¬ 
ing, Table IV and Eq. (E4) prove the general property 

0 for any two doublets A and A' in the CF spectrum of 
HTO. 


and (V'slil ^ are non zero for both states of the 

ground state doublet (m = 0,1). This is because in gen¬ 
eral the overlap between the perturbed A-states, in col¬ 
umn three and four, with the B-states, in columns five 
and six, is non-zero. Here, for brevity, only the results 
for ) are shown explicitly: 


b. B-doublets 

The other type of doublets in the CF spectrum of HTO 
are the H-doublets. The matrix elements V 
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= X! Jm( +a_(M-l)&-M 

M^2,5,8 


) = X! (-1)^Jm( «m&-(m-i) -a-(M-i)^M) ■ 
M^2,5,8 


(E5) 


Analougously one can show that also (V's^ I V \')pA^i) and 

(V's+il ^ IV'a+i) are non zero. 

Using their conjugation properties, one finds that 

^ I U = {i;f I U 

whose implications, in the context of the ground state 
splitting, are discussed in Sec. IIIB of the main text. 
Namely, the contribution to second order splitting due 
to i3-doublets vanishes identically. 

Whereas for notational convenience we have worked 
with a given choice of eigenstates for both the GS doublet 
and excited state doublets, the main results are indepen¬ 
dent of it. For instance, one can verify with a few lines 


of algebra that the relations in Eq. (E6) are invariant 
under generic basis transformations within each doublet 
involved. 


c. Singlets 

Another interesting feature of the CF eigenstates for 
HTO is the structure of the singlets and 

These are shown respectively in the fifth and sixth col¬ 
umn (from the left) of Table IV. To avoid confusion, it 
is important to underline that, in general, Si ^ s) for all 
i,j. The perturbative coupling of the singlets with the 
ground state doublet is non vanishing for both kind of 
singlets. Here, for brevity, we show explicitly only the 
matrix elements for 


U IV'®) = S3 (^a2j2 - a-4jl4j + S6 ya5j^ + a-Tjlyj -k soa_uli, 
IV'®) = 4 + a-4jl4^ + 4 “ a-rj-T^ ■ 


(E7) 


Analogously, it is straightforward to show that 

IV’i+i) 4 0 and I4a|i) 4 0. 

The matrix elements coupling the ground state dou¬ 
blet to the singlets provide the only non-vanishing 
second-order contribution to the ground state splitting in 
Eq. (11), marking the difference in the power law depen¬ 
dence found for HTO and DTO, as discussed in Sec. HI B. 

2. DTO 

All the energy levels in the crystal-field spectrum of 
DTO are doublets (see Fig. 2b). In Table V these are 
distinguished into A and B doublets, in analogy with 
HTO. 


a. A-doublets 

The two basis states, and of fh® ground 

doublet of the DTO crystal-field spectrum are of type A, 


as defined inJTable V. It is_then straightforward to verify 
that l4,o = 4,1 = 4.1 = 4,1 = 0 since in general 

u 1^4 f) = (iAi°Vil ^ I4i°|i) = (4i°4il V ) = 0, 

(E8) 

by comparing the first pair and the second pair of 
columns in Table V. The null matrix elements in Eq. (E8) 
are responsible for the vanishing of the first order contri¬ 
bution to the ground state splitting in Eq. (11). 


b. B-doublets 

As for HTO, also for DTO the matrix elements cou¬ 
pling the A and B doublets are non-vanishing: 
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15/2 

M=-9/2 
15/2 

(V'siil (-1)“’^® (jMaM&-(M-l) +J-Ma-MfcM-l ) , 

M=-9/2 


(E9) 


where the sum over the M quantum numbers runs, 
from —9/2 to 15/2, in intervals of 3 (M = 
-9/2, -3/2, 3/2, 9/2, 15/2). 

Their conjugation properties give: 

(<Wl<ii) = -(<Wl4°ii) ^ ^ 

(ElO) 

whose signs are opposite to the case of HTO in Eq. (E6). 
Similarly to the case of HTO however, Eqs. (ElO) give 
a vanishing second order contribution to the splitting of 
the ground state doublet. Since there are no singlets in 


DTO, no splitting at all takes place to second order. 

Because all the matrix elements in Eq. (E8) are null, 
the matrix elements in Eq. (ElO) are the only ones ulti¬ 
mately responsible for the DTO energy splitting, which 
takes place to third order in (transverse field) perturba¬ 
tion theory, as illustrated in Sec. III. 

We stress once again that, whereas for notational con¬ 
venience we have worked with a given choice of eigen¬ 
states for both the GS doublet and excited state doublets, 
the main results are independent of it. For instance, one 
can readily verify that the relations in Eq. (ElO) are in¬ 
variant under generic basis transformations within each 
doublet involved. 
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